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Abstract 

We study pion condensation at zero temperature in a hard-wall holographic model of hadrons 

O '. 

p^ , with isospin chemical potential. We find that the transition from the hadronic phase to the pion 

condensate phase is first order except in a certain limit of model parameters. Our analysis suggests 

f^ ' that immediately across the phase boundary the condensate acts as a stiff medium approaching 

the Zel'dovich limit of equal energy density and pressure. 



en 
rn 

O 

o 



X 



I. INTRODUCTION 



Dense nuclear matter generically carries net isospin and consequently has a nonvanishing 
isospin chemical potential, /i/. The Relativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) are presently the premier laboratories capable of studying the quark- 
gluon plasma phase of nuclear matter. The heavy nuclei in these experiments carry isospin, 
so the finite temperature system after collision has a nonvanishing isospin chemical potential. 
Neutron stars have still larger chemical potential, and the behavior of such objects depends 
on the phase structure of QCD in cold, dense environments ^miM- -^^ large chemical potential 
it is expected that mesons condense, beginning with the pions at /i/ ~ m^r (5H7|. Thephase 
diagram of QCD at high density has been explored via the chiral Lagrangian |8|, |9|, the 
Nambu-Jona-Lasinio model 1CH12|. and lattice QCD 13l-ll6|. The general consensus is 
that at low temperatures and vanishing baryon chemical potential there is a second order 
transition from the hadronic phase to a pion condensate phase at a critical isospin chemical 



potential /x/ around the pion mass m^ 



The AdS/CFT correspondence 



(holographic) models of QCD 
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40 MeV, or larger at finite temperature. 



17-[l9| has provided motivation for extra-dimensional 



28( 1 . Both explicit and spontaneous chiral symmetry 



breaking may be built into the extra-dimensional mode 



similar to extended hidden-local-symmetry models 29 



s, resulting in an effective description 



30|. The AdS/CFT correspondence 



maps sources and expectation values of field theory operators to backgrounds of extra- 
dimensional fields with corresponding quantum numbers. By studying fiuctuations about a 
prescribed background, the model makes predictions for field theory observables. 

We study the behavior of matter with isospin chemical potential in a holographic model 



of hadrons. Related analyses in other holographic systems appear elsewhere [e.g. Refs. 31 



34l|). We are motivated to study the hard- wall model of Refs. 23|, |24l in light of recent 



suggestions that the pion condensate phase is absent in that model [35|. In contrast, we 



identify the pion condensate 



The hard-wall models 
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Dhase and study its properties. 



24j are the simplest holographic models which capture certain 



features of QCD. The hard-wall geometry is a region of 4-1-1 dimensional (5D) Anti-de 
Sitter space preserving the isometrics of 3-1-1 dimensional (4D) Minkowski space. The main 
motivation for this choice of spacetime is its simplicity, although arguments have been made 
to support the choice of Anti-de Sitter space in light of asymptotic freedom at high energies 



and evidence for an approximate conformal invariance in QCD at lower energies |36l . |37] . 
The presence of a wall, which terminates the spacetime at what is referred to as the IR 
boundary, leads to a discrete spectrum of Kaluza-Klein modes identified via their quantum 
numbers with towers of hadronic resonances. 

Global symmetries are lifted to gauge invariances in the holographic description. In order 
to include the approximate chiral symmetry of the up and down quarks, SU(2)LxSU(2)i^ 
gauge fields are included in the 5D model. To mimic the pattern of chiral symmetry breaking, 
a set of scalar fields transforming in the bifundamental representation of the chiral symmetry 
group is introduced. The quantum numbers of the scalar fields are those of the scalar quark 
bilinear g^^gjj, where i and j are flavor indices which are now gauge indices. The scalar fields 
have a background profile which preserves the 4D Lorentz invariance of the spacetime but 
breaks the chiral symmetry to the diagonal isospin subgroup. 

The linearized equations of motion for the scalar field have two independent background 
solutions: a normalizable mode and a non-normalizable mode, the difference being that the 
normalizable mode has a finite effective 4D action, while the non-normalizable mode has a 
divergent action. If the non-normalizable mode is turned on, then the theory is modified so 
as to explicitly break the chiral symmetry, as if by a quark mass; if the normalizable mode 
has a nonvanishing background then there is a spontaneous breaking of the chiral symmetry, 
as if by a contribution to the chiral condensate (^^Q'^j). 

The chemical potential for the third component of isospin acts as a source for the isospin 
number density 

N3 = q\:Tlqi + ql^T^q^, (1.1) 

which is the time component of the isospin current 

Jr = tY^q', (1.2) 

where T" = cr"/2 are the generators of the isospin SU(2), and the subscript V represents 
the vector subgroup of the chiral symmetry. A source for the time component of the current 
couples as would the time component of a background gauge field. Hence, a non-normalizable 
background for the time component of the vector combination of 5D gauge fields mimics an 
isospin chemical potential in the 4D effective theory. 

As the magnitude of the isospin chemical potential is increased above the pion mass, 
it becomes energetically favorable for a pion condensate to form. We find that the phase 



transition is first order in the hard-wall model unless the 5D gauge coupling vanishes. The 
speed of sound Cs at high temperatures was conjectured to satisfy a "sound bound" c^ < 1/3 
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where c^ = 1/3 is the conformal limit. Fluctuations in the condensate at zero 
temperature violate the "sound bound," except near the phase boundary and then only if 
the 5D gauge coupling is small enough. Violation of the sound bound at low temperature 



is not unusual [39i] and has also been observed in certain D-brane systems 3^ and in a 



holographic model describing matter at a Lifshitz point 40|. 

To describe systems at nonvanishing temperature, extra-dimensional models are modified 
to include a black-hole horizon. However, we will focus on the zero-temperature phase of 
isospin matter, which corresponds to the original hard-wall background without a black- 
hole horizon. For simplicity we do not include chemical potentials except for isospin, so our 
analysis provides only a narrow cross section of the phase structure of the model. Extensions 
of these results to nonvanishing temperature and baryon chemical potential, and to include 



strange quarks and Kaon condensation 4l|], may shed light on the phases of matter in 



neutron stars and other extreme environments. 



II. HOLOGRAPHIC PION CONDENSATION 



The action for the 5D hard- wall model with chiral symmetry is given by 23|, |24| |. 

S = j <fx v/^Tr ||DX|^ + 3 |X|^ - ^ (F| + F^) | , (2.1) 

where DmX = OmX -HmX + iXRM, Lm = L^T" and Ff^^ = dMLN-duLN-i [Lm, Ln] 
(similarly for R), and we normalize the gauge kinetic term as in [23]. The spacetime in the 
hard- wall model is a slice of AdS^: 

ds = a{zY [ri^ydx^dx" — dz^^ , e < z < Zm, 

where a{z) = 1/z in units of the AdS curvature scale, and rj^^, is the 4D Minkowski metric 
with mostly negative signature. Greek indices range from to 3, and capital Latin indices 
from to 4, with x^ also denoted by z. The scalar fields X transform in the bifundamental 
representation of the SU(2)lxSU(2)r gauge invariance. 

Chiral symmetry breaking is provided by the background solution to the X field equation 
of motion, 

Xo{z) = -{m,z + az'') =^v, (2.2) 



where rriq is the quark mass matrix responsible for sourcing a, the chiral condensate. The 
bulk vector gauge field V^^ = l/2{L1,j + R1^) is dual to the isospin vector current operator. 
We work in the gauge L" = i?^ = 0. The linearized equations of motion for the transverse 



part of V^ are 



d. (;5.^;) - Id^d'^V^ = 0. (2.3) 



The background solutions for V^^ are of the form 

V,'{z) = c, + 'jz\ (2.4) 

where the coefficient of the non-normalizable mode, ci, is identified with the chemical po- 
tential for the third component of isospin ///; and C2 is proportional to the spontaneously 
generated background isospin number density, which we assume to vanish. Hence, the back- 
ground gauge field is uniform, 

V,' = fii. (2.5) 

The pions are identified with solutions to the linearized coupled equations of motion for 
the Goldstone modes in the scalar fields X, which mix with the longitudinal part of the 
axial vector field A'^ = {L'^^ — -R^)/2 = 9^0". We parametrize the Goldstone modes by fields 
tt" such that, 



(2.6) 



X = Xq exp [i2^°-T''] 

= Xq (cos b + i {n"'a"-) sin b) , 

where b = a/tt^tt^ and n^ = b^^n'^. The action (12. ip takes the form: 

S = f d^xy^hx^ (5a./(cos b)d^' (cos b) + Om (n"^ sin b) d^ {n" sin b) 

- 2/x/a-2(9o in" sin b) e^^^n" sin b - 2a~^d^ (n" sin b) cos b d^cj)" 

+ 2a-25^(cos b)n'' sin b d'^cj)" + 2^iia-^ cos b e"''^"rf sin b d^c^" 

+fi^ja-^ sin^ b n'n'^e'^'e'^^" + a'^ cos^ b d ^(j)" d^" (t)" + a'^ sin^ b 7fd^,<P''7fd^'<P") 

-«"'^ [/^7 (a,0^9>i + a,0^aV2) _ d^d.rd.d^r + O {{Alf)] |, (2.7) 

where contractions of Greek indices are done with 77^^, and those of capital Latin indices are 
done using the full metric Qmn- 



To quadratic order in tt" and 0", the action is 
The hnearized equations of motion for tt" and 0" are: 



/,o_^o 



„2, 






k^d, (a9,0O) 



mn0° - mn7r° = -^d. (v'^a^d^TT^) 



(2.9) 



m±0=^ - (m± =F /i/) 7r=^ = ^;^^<9^ (a^^^"^) , 



(m| =F fiim±) (p^ - (m| =F 2/i/m± + /if) 7r=^ = ^(9^ (t;2a3(9^7r=^) , 
where n^ = -^{n^ + 71^), n^ = ^(7r+ — vr^), and vr^ = 7r° (similarly for 0"). These equations 
are evaluated in the pion rest frame {q = 0), identifying the pion frequency with the effective 
pion mass in the isospin background. Working in the gauge A1^ = Aj^^ = 0, the fields satisfy 
the boundary conditions 9^0^'°(-2m) = 0^''^(e) = 7r^'°(e) = 0. The boundary condition at z^ 
corresponds to the gauge- invariant condition F^^{zm) = F^^{zm) = 0, but this choice is not 
unique and is made for simplicity. 

By eliminating 0^'*^ we obtain equations of motion for 7r^'° alone. 

d. (^9, {v'a^d^-K^)) + mld^-n^ - 9^0^8,71^ = 0, 
d. {^d, {v^a^d^Ti^)) + {ml T 2fiimi + /i?) 9,7r± (2.10) 

—glv'^a^dzTT^ = 0. 
Except for the replacement of the eigenvalue ttIq with (m^ =1= 2jj,jm± + jif^ , the fields tt^ 
and 7r° are solutions to the same differential equation with the same boundary conditions, 
with identical eigenvalues 



m, 



,2 _ /_2 ^ o. ™ 1.2 




miT^fiim^+i^j). (2.11) 



The neutral pion is unaffected by the chemical potential, so identifying ttiq with the pion 
mass in vacuum, m,r, we find a relation for the masses of the charged pions: 

m± = ±fj,j + rriTr, (2-12) 

where m,r = ""^o- For \f^i\ > f^-w a charged pion mass becomes negative, indicating the 
instability to pion condensation. 
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III. PROPERTIES OF THE PION CONDENSATE PHASE 

Since the pattern of chiral symmetry breaking is built into our holographic model, we 
expect to reproduce predictions of the chiral Lagrangian, at least qualitatively. The pion 
effective theory is determined by the action on the solution for the pion mode discussed in 
the previous section, integrated over the extra dimension. 

A. Decoupling the 5D gauge fields 

The limit (75 — )■ provides the most direct comparison to previous results. In that limit 
the fluctuations of the 5D gauge fields decouple from the pion physics. The corresponding 4D 
effective theory is similar to the chiral Lagrangian with isospin chemical potential included 
as a background for a 4D isospin gauge field, as in Ref. |8|. In terms of the unitary fields 



E = exp 
the leading order chiral Lagrangian is 



fn 



(3.1) 



£4^ = i^Tr (V.S VE^) + ^^Tr (S + S^) , (3.2) 

where VqS = SqS — z^ [cts, S] and Vj = (9,. Expanding to second order in the pion fields, 
the Lagrangian takes the form. 

The instability when |/i/| > m.,^ signals the phase transition to a pion condensate phase. 
Estimation of the value of the condensate and related observables requires an extension 
of the analysis to higher order in the pion fields, which we perform in the holographic 
description. 

By design, the analysis of the 5D model is similar to the chiral Lagrangian analysis above. 
In the limit g^ — )■ 0, we neglect couplings to the longitudinal gauge field (9^0'*. The action 
(12. 7p takes the form 

- 2fiia-^do {n' sin b) e"=^^n" sin b + fi^ja'^ sin^ b n'n'^e'^'e'^'^^ \ , (3.4) 



where b = a/tt^tt^ and n'^ = b ^tt*^ as before. The hnearized equations of motion for the pion 
fields are now, 



2_0,± 



1 



m^TT 



v'^a? 



d, (t;2a39,7r°'±) . (3.5) 



The condensate is a static configuration rotationally invariant in x^,x^,x^. The action on 
such configurations gives the condensate effective potential, 

K//,,3=o = jdz v{zfa{zf Q (^) ' + ^ sin^ b {^-^^ ' - ^ sin^ b n'^n' {5'^'' - 5'^H'^) ) . 

(3.6) 

The effective potential increases with \dn'^/dz\, so dn'^/dz = in the ground state. The profile 
of b{z) is determined from the solution to the equations of motion for the pion Kaluza-Klein 
mode. Expanding to fourth order in the pion fields, 

Veff^^^O = P dz V^Z' U (^) ' - /^? (-(-)'- ^ + ■■■ ) ^'^" {^'" - ^''n 

' dzv^z^ ^ (mln{zf - fijn'n'' {6''' - 6'^^6'^^) (n{zf - ^^ + ■■■))' (3-7) 

where we used the linearized equation of motion (13.51) in the last line. 

For l/i/l > rriTT it is energetically favorable to turn on the charged pions. The pion field is 
normalized by its kinetic term in the effective 4D theory, so we define 7c"'{z) = if{z) ti^ such 
that 

dzv'^a^^{zf = 1, (3.8) 

and vr" is the pion condensate (vr°). 

Minimizing Vejf expanded to (9((7r")^), we find that the transition is smooth (second 
order), and for ^ > m.,^ we obtain. 



where rj = J^"" dzv'^a?T^{zY. We then find. 



39/^ m. 



2\ 2 



K//,.5=o(vr^) = -^/^/ (^1 - ^ ) • (3.10) 



The isospin number density is 



dVeff 3/i/ / m. 



dfj,l 4:7] \ fj,j 



,4 



ni = -^^^=^-^ 1--^ . (3.11) 



We can express fj in terms of /^ in this model by the AdS/CFT determination of Jt^. 
The correlator of a product of axial currents has a pion pole at zero momentum transfer 
in the chiral limit, with residue equal to /^. The AdS/CFT correspondence determines the 
correlation function in terms of a bulk-to-boundary propagator which solves the linearized 
equations of motion for the transverse part of the axial vector field. For more details in the 



context of the present model, see Refs. [23|, |2J]. We summarize the results here. 



The linearized equation of motion for the transverse part of the axial vector field Al^{q, z) 
is given by, 

d,{ad,A;) + ^-A';,-v'a'glA';] =0. (3.12) 

The bulk-to-boundary propagator A{q, z) describes the solution to (13.121) of the form 
y4^(g, z) = A{q,z)AQ^{q), with boundary conditions dzA{q,z)\z^ = and A{q,e) = 1. In 
terms of the bulk-to-boundary propagator the AdS/CFT prediction for the pion decay con- 
stant is 



^^_l_dzA{0,z) 



(3.13) 

95 z 

If (75 = then the bulk-to-boundary propagator at g^ = is uniform, A{0,z) = 1. To 
next order in g^, we obtain 

-dzA{0, z) = -gl / dzv{zr/~z' + 0{gl/zt). (3.14) 

From (I3.13P we obtain in the (75 — )■ limit. 



p^ = j:-dzv{z)yz' 

rriqazl^ + ml \og{zm/e). 






(3.15) 



Note that in the absence of boundary counterterms we must choose e such that log{zm/^) <^ 
az'^/niq to respect the chiral limit. We choose e = 1/(10^ MeV). In that case the integral 
defining /^ is dominated by the region where the pion wavefunction Tf{z) is approximately 
constant. Comparing (I3.15P with (13.81) . we learn that the pion wavefunction tt{z) ^ 1//^ 
except for a region of small z, as in Figure [H Similarly, from the integral defining fj we have 

V ^ l/fl (3.16) 

For a concrete example, fixing the mass of the lightest KK mode of the vector field V^ to 776 
MeV determines Zm = 1/(323 MeV) [23]. Then with mq=4.25 MeV and chiral condensate 




FIG. 1: Pion eigenfunction with niq = 4.25 MeV, a = (263 MeV)^, and Zm = 1/(323 Mev). 

a = (263 MeV)3 we find physical values m^=140 MeV and A=92 MeV in the 5-5-^0 limit. 
With these values of the parameters, we find 

dzv{zya{zfn{zy = 1/(91 MeV)^, (3.17) 

which is approximately 1//^ as expected. Note that the Gell-Mann-Oakes-Renner relation 
is approximately satisfied, m^/|/(2m^o") = 1.07 ~ 1. 

We now have the holographic prediction of the equation of state: 



For comparison, the corresponding prediction from the 4D chiral Lagrangian (13. 2 p is [S' 



(3.18) 



niD = /^/i/ (1 f / ' 



m: 



(3.19) 



which differs from the holographic prediction by an overall factor of 4/3. This overall 
factor drops out of the ratio of pressure to energy density and the speed of sound at zero 
temperature. The number densities are plotted in Figure [2j The model is not expected to 
be valid for /i/ > rrip ^ 5.5771^^, but we plot the model prediction here and below over the 
entire range of fij. 

The pressure p and energy density e in the pion condensate medium are determined by 

ii) urn- 



niini 



Pif^i 



MI 



rij djl 



8/i? 



nj 



^M = I f^i dh = -^ {fij - ml) {fij + 3ml) 
'0 o/^/ 



(3.20) 
(3.21) 
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FIG. 2: Isospin number density. The bottom red curve is the prediction of the hard-wah model 
with rug = 4.25 MeV, a = (263 MeV)^ and Zm = 1/(323 MeV). The top blue curve is the result 
from Ref. ^ quoted in (l3T9]l . 



This gives 



and 



P 



li]-ml 



e ix] + 3ml ' 



2 dp nj- m% 



(3.22) 



(3.23) 



' de fij + Smf 
The speed of sound violates the sound bound c^ < 1/3 except near the phase transition 

boundary at fij = 711^^. 



B. Pion condensation with dynamical 5D gauge bosons 

Having understood how the fifs — !■ limit of holographic QCD reproduces the qualitative 
behavior of isospin matter at low temperature expected from the chiral Lagrangian, we 
now consider the more general situation including couplings to the 5D gauge fields. In the 
calculations below we will take g^ = 2-7?, which makes the holographic prediction of the 
vector current polarization at large momentum transfer agree with perturbative three-color 
QCD 



23 



24| . 



We first construct an approximate solution to (12. 9p as an expansion in m.,^, as in Ref. 23 1 
Combining the first two equations of (12. 9p we get 



m^dz(, 



v'^a^gld^'K^. 



(3.24) 
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Recalling the boundary conditions 7r(e) = 0(e) = dz(t)\zm = 0; to zerotli order in m.„ the 
solution is 7r° = 0, and 0° satisfies the same equation as the bulk-to-boundary propagator 
in fl3.12p . so we set 0°(;z) = A{0,z) — 1. Away from the boundary at 2; = e consistency of 
the approximate solution with the first of the equations in fl2.9l) requires 7r° ^ —1. Recalling 
that the charged pions have the same wavefunction as the neutral pion, we temporarily 
normalize the fields so that tt^{z) = tt^{z). Then 0^ = (1 =]= ^i/m±) [A(0, 2) — 1]. 
Integrating the action (12.81) by parts we get, 



S 



d xdz 



v'^a^ (tt+tt" 



vr 



vr 



+ 0+0 ) 2^2(0^2 



d^7r^{x)d''7r-{x) 



+ [dz (v'^a^dzTv^) TV + jj'jv'^a^'K^'K ] 7r+(a;)7r (x) 



vr X vr [X 



(3.25) 

where functions without an argument are understood to be functions of z. Ignoring the vr° 
terms, which can be obtained by taking /i/ — )■ 0, we can use the third of Eqs. (12. 9p to solve 
for 0+ and obtain 



S 



d X 



/^? 



+ 






1 



dzv'^a^Tm ^ / dzdyiady. 



IT 



m^gl 



dz dz{adz 



vr 



9^vr+(x)9^vr-(x) 



dzdz{adz(l)'^)TT + fij I dzv'^a'^mi 
1 



dz f a vrvr + 



9l 



~i 






(iz t> a vrvr + 



dzdz{adz<t> tt) 
1 



vr+(a;)vr (x) 
9tvr^(x)7r+(x) 



^i 



dz dz{adz(j) )tt 



9ivr+(x)vr-(x) L (3.26) 



where we have used vr = vr'^(2) = tt^{z). We now use the approximate solutions for 0+ and 
make the approximation vr = — 1 in those integrals dominated by the region where the pion 
wavefunction is fiat. Defining a = -k I dz v'^a? vrvr and making use of equation (I3.13p . we get 



S 



d X 



ml — ajji] 
ml- fij 



f%n+{x)d'^n~{x) - K - «-"?) /V(a:)vr-(x) 



-2zf,jf( '^\ ''i )d,n^x)n-{x)]. (3.27) 
The resulting Lagrangian after canonically normalizing the kinetic term is, transforming 



12 



back to (7r\7r^,7r^), 

This agrees with the leading order 4D chiral Lagrangian (13.31) . 

The effective potential for static configurations of 7r(x) takes the same form as ( 13. 6p 
because the longitudinal gauge bosons are derivatively coupled. We make the ansatz that 
rr" = 0; the pion expectation value is only in the {tt^,tt'^) plane. With the approximate 
solutions for tt{z) and (f){z) described above, and applying the canonical rescaling of the 
kinetic term, the effective potential at 0{b'^) becomes 

V.„ W-lhl- .5) ^^ + i.5/- (£ff|j) \^\ (3.29) 

where rj = j2 J dzv'^o? tttttttt, a = -4- J dzv'^o? nn, and /^ is given by fl3.13p . Note that tt{z) 

J TV J IT 

here is the approximate solution for small m.^- described earlier, and is not the canonically 
normalized field. We can now solve for the minimum of the potential, dVeff/db = 0: 



3 P. K - c^f^lf 
2r]nj (/i?-m2) 



bl = ^^ "'7^2 Zii ■ (3-30) 

The isospin number density is 



dVeff „9 1 3 / 9 m. 



,4 



ni = —j^ = P^ift-l «^ - ^ • (3.31) 

The phase transition is first order, with the order parameter jumping at the phase boundary 
/i/ = 171.,^. The minimum value of the free energy Ve// at (9(6^) is discontinuous at the critical 
point /i/ = m,r5 but the perturbative expansion breaks down near the critical point and we 
expect a nonperturbative analysis to confirm continuity in the free energy across the phase 
boundary. However, we can say with confidence that the transition is not second order in 
this model, because if it were then the order parameter b would vary smoothly and we would 
expect a perturbative analysis to be valid. Although perturbation theory breaks down near 
the transition, we plot the perturbative prediction for the number density in Figure |3l 
The ratio of the pressure to energy density is now 

p ^ a'^fij - ml 
e Q^ii] + 3m2 ' 
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FIG. 3: Isospin number density. The top red curve is the perturbative prediction of the hard- 
wall model with parameters fit to rup, m-j^ and f-j^: nig = 2.26 MeV, a = (333 MeV)'^, and 



Zm = 1/(323 MeV). This set of parameters gives a = 3.66 and ij = 3.60, in (j3.3ip . The bottom 
blue curve is the result from Ref. 8| given in (j3.19p . 



and 



2 4 4 



This is plotted next to the chiral Lagrangian prediction in Figure HI Note that the speed 
of sound exceeds the sound bound c^ = 1/3 throughout the pion condensate phase at zero 
temperature. 



IV. A COMMENT ON THE CHIRAL LAGRANGIAN 



Leading order chiral perturbation theory predicts that the transition to the pion con- 
densate phase is second order. We have learned that gauging the chiral symmetry in the 
holographic model qualitatively modifies predictions for pion condensation at zero tempera- 
ture. The transition becomes first order, and the medium becomes stiff immediately beyond 
the phase boundary. Including higher derivative terms in the chiral Lagrangian can have 
similar consequences, as we will now demonstrate. Consider the Lagrangian 

C = f Tr [D^SD^St] + Ili/^Tr [S + St] + ai (Tr [D^.SD'^St])^ 
+a2Tr [D^SLi^St] Tr [D^SLi'^St] , 



(4.1) 
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m„ 200 300 400 500 600 700 "V 



y"/ (MeV) 



FIG. 4: Speed of Sound. The upper red curve is the perturbative prediction of the hard- wall model 
for the speed of sound with ruq = 2.26 MeV, a = (333 MeV)^, and Zm = 1/(323 MeV). This set 
of parameters gives the value a = 3.66 in (j3.33p . The bottom blue curve is the prediction based 
on (|3.19p . The top and bottom dashed lines represent the speed of light and the conformal limit 
c2 = 1/3. 

where for this analysis ai and 02 are arbitrary parameters. Once again we take the static 
part of the Lagrangian to get an expression for the effective potential. Defining S = cos b + 
i {n"'a"') sin 6, we have 



Ve ff {cos b) 



(l-cos^fe) (1 



n^n^ 



)-mlJtcosb-a,^{l 



cos 



6)^(1 



n n 



r. 

(4.2) 



where Oi = ■S{c(i+a2)- At the minimum of K//, n^ = 0, and we find a region of ai parameter 

J IT 

space where the phase transition is first order. That is, as ai increases past a critical value 
^crit _ 1/(2777,2)^ the phase transition changes from second to first order. This is illustrated 
in Figure [5] . However, f^af^^ = 0.22 is much larger than the typical low energy coefficients 
in the chiral Lagrangian inferred by experiment (/i(m7r) = (—4 ± 6) x 10^^, hiiTi,^) = 



(9.1 ±0.2) X 10' 
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43|. 



V. CONCLUSIONS 

We have studied pion condensation at zero temperature and finite isospin chemical poten- 
tial in a hard-wall model of holographic QCD with chiral symmetry breaking and massive 
pions. At the critical point fij = m^ the pion condenses, and our perturbative analysis 
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0.0 



0.2 



0.4 0.6 

cos(b) 



0.8 



1.0 




0.4 0.6 

cos(b} 




(b) ai critical. 
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cos(b) 
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FIG. 5: Each plot shows the phase transition for a different value of the oi parameter. The critical 
value of fij for pion condensation depends on oi. The three curves shown in each plot correspond 
to fij < fic, fJ'i = fJ-c, and /i/ > fic (top, middle, and bottom curves, respectively). Plot (a) shows 
the transition for ai less than the critical value. Plot (b) is with ai the critical value, while ai of 
(c) is larger. These plots assumed m.,^ = 139 MeV and fj^ = 92.4 MeV. 

suggests that the condensate creates a stiff medium approaching the Zel'dovich equation of 
state p = e. Sound propagation exceeds the conformal sound bound c^ = 1/3, except near 
the phase transition boundary if the 5D gauge couphng is small enough. The low-energy 
effective theory for pions as derived from the hard-wall model indicates that the transition 
from the hadronic phase to the condensate phase is first order, except in the limit of van- 
ishing 5D gauge coupling. This is in contrast to leading order chiral perturbation theory, 
which predicts a second order transition SL and lattice simulations which also seem to be 
consistent with a second order transition [14]. We have shown that even in chiral pertur- 
bation theory the transition can become first order if higher derivative terms in the chiral 
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Lagrangian have large enough coefficients. 
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